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Protograph-Based LDPC Code Design 
for Shaped Bit-Metric Decoding 

Fabian Steiner, Georg Bocherer and Gianluigi Liva 


Abstract —A protograph-based low-density parity-check 
(LDPC) code design technique for bandwidth-efficient coded 
modulation is presented. The approach jointly optimizes the 
LDPC code node degrees and the mapping of the coded bits 
to the bit-interleaved coded modulation (BICM) bit-channels. 
For BICM with uniform input and for BICM with probabilistic 
shaping, binary-input symmetric-output surrogate channels 
for the code design are used. The constructed codes for 
uniform inputs perform as good as the multi-edge type codes of 
Zhang and Kschischang (2013). For 8-ASK and 64-ASK with 
probabilistic shaping, codes of rates 2/3 and 5/6 with blocklength 
64800 are designed, which operate within 0.63 dB and 0.69 dB of 
| log 2 (l + SNR) for a target frame error rate of 10 -3 at spectral 
efficiencies of 1.38 and 4.25 bits/channel use, respectively. 


I. Introduction 

Bit-interleaved coded modulation (BICM) combines high 
order modulation with binary error correcting codes [1], [2], 
This makes BICM attractive for practical application and 
BICM is widely used in standards, e.g., in DVB-T2/S2/C2. 
At a BICM receiver, bit-metric decoding (BMD) is used [3, 
Sec. II], Achievable rates for BMD were investigated for 
uniformly distributed inputs in [3] and for non-uniformly 
distributed bits in [4], These results were generalized to non- 
uniformly distributed input symbols in [5] and [6], BICM 
with probabilistic shaping can achieve rates very close to the 
capacity l log 2 (l + SNR) of the additive white Gaussian noise 
(AWGN) channel [6]. The aim of this work is the design 
of low-density parity-check (LDPC) codes for BICM with 
probabilistic shaping to approach ilog 2 (l + SNR) at high 
spectral efficiencies. 

The key aspect of designing LDPC codes for BICM are the 
unequal error protection of the LDPC coded bits and the BICM 
bit-channels that are different for different bit-levels. A first 
approach is to take an existing LDPC code and to optimize 
the mapping of the coded bits to the BICM bit-levels. This 
was done, e.g., in [7]—[9], A more fundamental approach is 
to directly incorporate the different bit-channels in the code 
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design. This is done in [10], where the authors use multi¬ 
edge type (MET) codes [II] to parameterize the different bit- 
channels. They then extend the extrinsic information transfer 
charts (EXIT) [12] to multiple dimensions to design codes 
for quadrature amplitude modulation with 16 signal points 
(16-QAM). As 16-QAM can be constructed as the Cartesian 
product of two four point amplitude-shift keying (4-ASK) 
constellations, two different bit-channels are apparent. For 
constellations with more than two different bit-channels, the 
authors of [10] observe long runtimes of their multidimen¬ 
sional EXIT approach. Therefore, they suggest a high-order 
extension based on nesting, i.e., starting from m = 2 bit-levels, 
they successively extend their codes from m to m +1 bit-levels 
by optimizing in each step only the additional bit-level. 

In this work, we follow [10] and jointly optimize the code 
structure and the mapping of the coded bits to the BICM 
bit-levels, but propose a protograph-based design [13]. In 
[14], protographs are designed for coded modulation by a 
variable degree matched mapping (VDMM). Each bit-level is 
associated to a specific protograph variable node following 
the waterfilling approach (i.e., assigning the most protected 
coded bits to the bit-levels with highest bit-channel capaci¬ 
ties). Recently, a protograph-based coded modulation scheme 
was introduced in [15] by performing a one-to-one mapping 
between the constellation symbols and the codeword symbols 
of a non-binary protograph LDPC code. This requires the 
constellation order to match the field order on which the LDPC 
code is constructed. To the best of our knowledge, none of the 
above-mentioned approaches leads to a joint optimization of 
the binary LDPC protograph ensemble and the bit-mapping. 

We propose the following code design. We represent each 
bit-channel by a biAWGN surrogate. Our surrogates reflect 
both the BICM bit-channels and the input distribution. For 
the surrogate channels, we employ protograph EXIT (P-EXIT) 
analysis [16] to determine the ensemble iterative convergence 
threshold by accounting for the different bit-channels asso¬ 
ciated with the protograph variable nodes. The protograph 
ensemble is optimized with respect to the threshold by dif¬ 
ferential evolution. We optimize codes both for uniformly 
distributed inputs and for the probabilistic shaping scheme 
proposed in [5], Our optimized codes perform as good as the 
codes by Zhang and Kschischang [10]. Moreover, our design 
approach can be applied to arbitrary large constellations with 
shaped inputs. Our codes of rates 2/3 and 5/6 with blocklength 
64800 operate within 0.63 dB and 0.69 dB of ^ log 2 (l + SNR) 
for a target frame error rate of 10 3 at spectral efficiencies of 
1.38 and 4.25 bits/channel use, respectively. 

This paper is organized as follows. In Sec. II, we review 
LDPC coded modulation with probabilistic shaping. Sec. Ill in- 
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Figure 1. 8-ASK constellation with BRGC [17] labeling. 
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where H(-) denotes entropy. 

Remark 1. If the bits Bi, B 2 ,..., B m are independent, then 
(5) can be written as [3], [4] 

m 

7 ?bmd = ^ I (Bj; L,) ( 6 ) 

i=l 

where !(•; •) denotes mutual information. 


B m -► ( PL m |B m )-*- L m 

Figure 2. m parallel bit-channels that are different for different bit-levels. 

traduces the concept of protographs and explains the necessary 
adaptions for P-EXIT analysis. We present our code design 
approach in Sec. IV. In Sec. V, we discuss the performance 
of our codes for uniform and shaped inputs, respectively. We 
conclude in Sec. VI. 

II. Preliminaries 

A. Channel Model 

Consider the discrete time AWGN channel 

Y = AX + Z (1) 

where the noise term Z is zero mean, unit variance Gaussian, 
and where the input X is distributed on the normalized 2 m - 
ASK constellation 

= {±1, ±2,.... ±(2 m — 1)}. (2) 


B. LDPC Coded Modulation 

1) LDPC codes: LDPC codes are linear block codes with a 
sparse (n— k) x n parity-check matrix H. The matrix H can 
be represented by a Tanner graph [18, Section 3.3] consisting 
of variable nodes V t , i £ {1,... , 71 } and check nodes Cj,j £ 
{1,... ,n — k}. In the following, the set Af(Vf) denotes all 
neighbors of variable node Vi, i.e., all check nodes Cj with 
a one occuring in the i-th column of H ; the set Af(Cj) is 
defined as all variable nodes V] possessing a one in the j- th 
row of H. The variable and check node degrees of the *-th 
variable and j-th check node are then given as the cardinality 
of the sets A f{Vf) and Af(Cj), respectively. 

2 ) Transmitter: The canonical model of a BICM transmitter 
[19, Sec. 2.3] is displayed in Fig. 3a. k bits are encoded by 
a rate k/n encoder. The coded bits are then permuted by an 
interleaver and a mapper modulates the bits to a sequence of 
channel input symbols of length rtc = n/m. We adapt this 
general BICM model as follows. We describe the mapping 
function by specifying to which bit-level each coded bit is 
mapped. The mapper input is 

(Ba-.-B^B? 


The constellation spacing A controls the average power 
E[|AX| 2 ], where E[-] denotes expectation. The signal-to-noise 
ratio (SNR) is SNR = E[|AX| 2 ]/l. Each signal point x £ X 
is labeled by to bits B = (Bi, B 2 ,..., B m ), where B; denotes 
the i-th bit-level. Throughout this work, we label by the binary 
reflected Gray code (BRGC) [17], e.g., see Fig. 1. Let py|B be 
the transition function of the memoryless channel with input 
B and output Y. At the receiver, a soft-demapper calculates 
for each bit-level i the soft-information 
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channel likelihood 



input priors 


(3) 


which can be calculated by means of the joint distribution 


PBiib^PYiBiiylh) = PY\B{y\ a ) p B{a). (4) 

aG{0,l} 771: cn=bi 


A bit-metric decoder uses the soft-information L\, L 2 , . ■ ■, L m 
to estimate the transmitted data. For the decoder, the channel 
appears as m parallel bit-channels, see Fig. 2. Bit-metric 
decoding can achieve the rate [6, Theorem 1] 


i?BMD = H(B)-^H(B ; |Y) 
1=1 


m 

H(B)-^H(B l |L i ) 


i =1 

(5) 


where 

(B2 • • • B m )^ = (B2 • • • B m )i(B 2 • • • B m ) 2 • • • (B2 • • • B,,,)^ 

and 

B? = B n B 12 ---B lrk . 

We put the bits for bit-level 1 in the end. At this point, this 
choice is arbitrary, but it will become useful later on when we 
consider shaped inputs. The memoryless modulation function 
of the mapper is 

BuB 2i ■ ■ ■ B mi i-A XB li B 2i -B mi ■ (7) 

For LDPC coded BICM, each coded bit has a variable node 
degree and it is transmitted over one of m bit-channels. 
The resulting mapping of variable node degrees to bit-levels 
strongly influences the error performance of the system [7], 
If the LDPC code is given and fixed, we have to choose the 
interleaver appropriately. Since we are going to design new 
LDPC codes, we instead remove the interleaver and account 
for the mapping of variable node degrees to bit-levels in our 
code design. We choose a systematic rate (m — 1 )/m LDPC 
encoder, which realizes the mapping 

(B 2 • • • Bm)' ^ (Ba-'-B^B?. (8) 

Our model can be extended to code rates higher than (m — 

1 )/m, see [21]. For clarity of exposure, we only consider rate 
(to — 1 )/to codes in the present work. 
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(a) Canonical BICM model [19, Sec. 2.3]. 

Matcher j —> (B 2 ■ ■ ■ Bj 71 ) 77c — LDPC Encoder j —* (B 2 • ■ ■ B m ) ?7c Bj iL — Bit-Mapper j —> X 7 ^ 

(b) BICM with probabilistic shaping [5], 

Figure 3. Fig. (a) shows the canonical BICM model of [19, Sec. 2.3], whereas Fig. (b) depicts the adapted model needed for probabilistic shaping as proposed in 
[5]: Independent uniformly distributed data bits LL = U 1 U 2 ■■ ■ U/, ;| are matched to n c strings (B 2 B 3 ■ • • B m );, i = 1 . 2 ..... n c that are distributed according 
to • The systematic LDPC encoder appends ric check bits = B 11 B 12 ■ ■ ■ B] r , c . The bit mapper maps this bit stream to signal points according to 

the BRGC labeling. The overall rate is L/ry [bits/channel use] and the rate of the LDPC code is (m— 1 )n c /(rn/n L ) = (m— 1 )/m. At the receiver, a bit-metric 
decoder calculates an estimate of (B 2 ■ ■ • B m )" L B 7t . A data estimate U /,!l = U 1 U 2 • • • U/, ;i is obtained by passing the estimate of (B 2 • • • B„,)'" c through a 
dematcher. For the matcher and the dematcher, we use [20], 


3) Shaped Inputs: We use the shaping scheme proposed in 
[21], The capacity-achieving distribution of ASK constellations 
in Gaussian noise is symmetric around the origin. Conse¬ 
quently, it induces a distribution / j b, b 2 ---b„, on the BRGC 
labeling with the following properties: 

• Bit-level Bi decides on the sign of the transmitted 
constellation point, see Fig. 1. Because of symmetry, 
bit-level Bi is uniformly distributed. 

• Bit-levels (B 2 • • ■ B m ) and bit-level Bi are independent. 
Bit-levels B 2 , ■ . ., B rn are correlated. 

The scheme [21] mimics the capacity-achieving distribution 
in the following way: first, generate bit-levels B2 • • • B m ac¬ 
cording to Pb 2 - B m , e.g., by using a distribution matcher 
(see Fig. 3b). In our implementation, we use the distribution 
matcher proposed in [20]. The systematic encoder copies the 
bits (Bo • ■ ■ Bjn,)^ to its output and leaves their distribution 
un-changed. The encoder appends check bits B” 7 that are 
approximately uniformly distributed because each check bit 
is a modulo two sum of many information bits [22, Sec. 7.1]. 
The signal point XB 1 B 2 - B rn selected by the bit-mapper then 
has approximately the capacity-achieving distribution. 

III. Protographs and EXIT Analysis for Parallel 
Channels 

A. Protograph-Based LDPC codes 

LDPC code ensembles as introduced in II-B1 are usually 
characterized by the degree profiles of the variable and check 
nodes. For instance, X(x) = J2dL 1 AdX d ~ 1 and p(x) = 
Yhd=iPd xd l are the edge-perspective variable and check 
node degree polynomials with maximum degree d v and d c , 
respectively. However, the degree profiles do not allow to 
characterize the mapping of variable nodes to the m different 
bit-channels resulting from our adapted BICM transmission 
scheme. In the following, we use protographs to incorporate 
the bit-mapping in our threshold analysis. 

Parity-check matrices are constructed from protographs as 
follows. Starting from a small bipartite graph represented via 
its basematrix A = [aik] of size M x N, one applies a copy- 
and-permute operation (also known as lifting) to create Q 
instances of the small graph and then permutes the edges so 
that the local edge connectivity remains the same. The Q repli¬ 
cas of variable node 14,fc £ {1 ,...,W} must be connected 
only to replicas of the neighbors of 14 while maintaining the 


Ci C 2 



C, C 2 C[ C’ 2 



Figure 4. Above, the Tanner graph of basematrix A with M = 2 and N = 3 
is displayed. Below, an example lifting with Q = 2 instances of the protograph 
is shown. 


original degrees for that specific edge. The resulting bipartite 
graph representing the final parity-check matrix H possesses 
n = Q ■ N variable nodes and n — k = Q ■ M check nodes. 
Parallel edges are allowed, but must be resolved during the 
copy-and-permute procedure. An example protograph with the 
corresponding basematrix and an example lifting for Q = 2 
are shown in Figure 4. 


B. Protograph EXIT 

The performance of protograph-based LDPC codes over 
parallel channels can be easily analyzed in the infinite block- 
length limit by suitably modifying the P-EXIT. We consider 
next the case where the codeword bits corresponding to the 
N protograph variable nodes are transmitted over N binary- 
input symmetric-output parallel channels. In the following, 
we denote by I y k ->Ci t ^ le mut ual information (MI) between 
the message sent at iteration l by the fc-th variable node 
to the /-th check node and the corresponding codeword bit. 
Similarly, I^ y denotes the MI between the message sent 
at iteration £ by the /-th check node to the fc-th variable node 
and the corresponding codeword bit. We further express the MI 
between the k -th channel output and input as l( h . The evolution 
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( 10 ) 


( 11 ) 


of the MI can be tracked by applying the recursion 


|C,i _ rV 
Wk^Ci ~ T fc 


E,j 

Vu 

EJ. 

Ci^V k 


= fi 


' ict 

,i y'c^-Vk’ 'k 

u (if 


v-4-a) 


with 


and 


|E.< 

C —VI 


■V->Ci 


:= (%. 
■- Os 


|E.< 

Ci—>Vfe > 'Ca-s-V* 




C M ^Vk , 


t E,e ,e,i 

Vi->Ci 1 Va-J-Ci > ■ ■ • > 'Viv-i-Ci , 


( 12 ) 

(13) 


where by convention we set l E ^ c = 1^! = 0 if 

Ci N (14). In (12), (13) we introduced the variable and 
check extrinsic information transfer (EXIT) functions f kl and 
f E ; , whose form depends on the underlying channel model. 
We moreover denote by l^ PP '‘ the MI between the a posteriori 
probability (APP) log-likelihood ratio (LLR) message com¬ 
puted at the fc-th variable node in the f-th iteration, and the 
corresponding codeword bit. Note that l^ PPP is a function of 

iM -1 nnf ) |ch 

'e-tv* and fe • 

Following [23], we define the protograph convergence re¬ 
gion C\ as the set of channel MI vectors l ch = (if 1 , if 1 ,..., if 1 ) 
for which l^ PP < converges to 1 for all k £ [l,iV], as 0 —> oo, 


The evaluation of (12) and (13) at the 0.-th iteration is given 
by (9) and (10), with |^ pp ^ evaluated as in (11) and where 
J (<7 c h) is 


J (cr c h) := 1 - 


(—a./Q 2 


V a 


7 T(T 


log 2 (l + e *) dz 


ch 


Example 2 (Convergence region of parallel binary erasure 
channels). If the underlying parallel channels are binary 
erasure channels (BECs) with erasure probabilities e*,, k 
1 .Y. the convergence region in terms of erasure proba¬ 

bility is denoted by 


C f := 


eii£2, ■ 


i ejv 


lfc PP ’ / -> 1, Vfc, 0 


Here, the recursion (12), (13) at the 0-th iteration is given by 

M 


E,i 

W-s-C, 


= 1 - 


A iE,£—1 TT ft \ Qi ' 

ek i 1 'cwwj ii^ 1 'cwvfcj 


/'=1 
I‘¥=i 


E,e 

Ci^Vk 


fvur' ft fvCcr 


k ’=1 

fcVfc 


C,ch := 


|l ch l£ PP ’^ —> 1, Vfc, 0 —> oo|. 


We will sometimes restate the convergence region with re¬ 
spect to other parameters characterizing the channel nuisance. 
Hence, we will denote the convergence region with respect to 
the generic channel parameter £ as Cf. 

Example 1 (Convergence region of parallel binary-input addi¬ 
tive white Gaussian noise channels). The convergence region 
of N parallel binary input additive white Gaussian noise 
(biAWGN) channels, with noise variances crj? h k and inputs 

alphabets j±of h fc /2|, k = 1,..., N, may be expressed with 
respect to the noise standard deviation as 


a ch := 


<Tch,l> ^ch,2? • • • 5 &ch,N 


APP,^ 


1 , Vfc, £ —> oo 


with 


M 


,APP,£ _ 

] k 


i ’=t 


E,e-i 




IV. Code Design via Surrogate Channels 

The m bit-channels in Fig. 2 are not output-symmetric 
and the input is not uniformly distributed when we perform 
shaping. Optimizing LDPC codes directly for the original bit- 
channels in Fig. 2 is therefore difficult. To overcome this 
problem, we use the paradigm of universal codes. A code is 
universal, if it performs similar on different channels, under 
certain conditions. Universality allows to replace the original 
channel by a surrogate channel for which code design is 
easier. We first derive a universality criterion and then construct 
surrogate channels. 
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Figure 5. Search trajectory for protograph optimization for uniform and shaped 
inputs. 


A. A Universality Criterion 

We start by considering the rate backoff criterion [21, 
Sec. VIII.C]. Let R* be the rate that can asymptotically be 
achieved by the considered code ensemble on the original 
channel and let R be the actual transmission rate. The rate 
backoff is 

R* - R. (14) 

In [21, Sec. VIII.C], it is assumed that the error performance 
is a function of R* — R. Consider our transmission scheme of 
Fig. 3. Suppose the code rate is c so that (1 — c)m bits per 
channel use consist of uniformly distributed redundancy bits 
on average and suppose further that the input distribution is 
Pb- Then, by (5) and Sec. II-B, the achievable rate and the 
transmission rate, respectively, are given by 

m 

R* =H(B)-^H(B i |L i ) (15) 

i —1 

R = H(B) - (1 - c)m. (16) 

Taking the difference, we have 

m 

iJ*-i2=(l-c)m-53H(B i |L i ). (17) 

i=1 

Note that the term (1 — c)m does not depend on the statistics, 
so that the sum of the conditional entropies determines the 
rate backoff. Since we want to account for the different bit- 
channels in our code design, we strengthen the rate backoff 
criterion and use as criterion the bit uncertainties 

U = { H(B i |L,)}™ 1 . (18) 


B. Surrogate Channels 

Suppose we optimized a code on some surrogate channel. 
If the performance of the code on the original channel cannot 
be improved by changing the bit-mapping, we succeeded to 
jointly optimize the code and the bit-mapping. In this case, 
we say that the surrogate channel is proper. If in contrary the 
performance of the optimized code on the original channel 
can be improved by changing the bit-mapping, our joint 
optimization attempt failed. We will use the following two 
surrogate channels for code optimization. 

1) BEC Surrogate: We replace the original transmission 
scheme by a scheme where all input bits B, are uniformly 
distributed and where each bit-channel P^_ ,g is a binary 
erasure channel (BEC) with erasure probability e*. Criterion 
(18) requires 

=e i =H(B i |L i ), i = l,2,...,m. (19) 

2) biAWGN Surrogate: All input bits B, are uniformly 
distributed and transmitted over biAWGN channels 

U = x§. + Z i (20) 

where Xg £ {±cr/ h i/2} and the noise Z, has variance er/ h ... 
Applying criterion (18) yields 

H(Bi|Li) =H(B,|Li), t = l,2,...,m. (21) 

This can be achieved by numerically solving for <J C b t i. 

C. Iterative Decoding Threshhold 

Since irregular LDPC codes have shown better performance 
than regular LDPC codes [24], we want to allow different 
variable node degrees for each of the m distinct bit-channels 
as well. In order to have up to D different degrees per bit- 
channel, the protograph matrix A must have at least N = D m 
variable nodes. We introduce a mapping function of the form 
T(k) = \k/D] to relate each variable node 14 with a 
corresponding bit-level T(k) £ {1,, m}. We next define the 
iterative decoding threshholds for uniform and shaped inputs, 
respectively. 

1) Uniform Inputs: For uniform inputs, the bit uncertainty 
set 

U = = {H(B,|Ax b + Z)}i=i (22) 

is a function of the SNR via the constellation scaling A. Let 
£W(SNR) denote th e vector of surrogate channel parameters as¬ 
sociated with the uncertainty set W(SNR). Given a protograph 
ensemble, its iterative decoding threshold is given by 

minSNR s.t. £ W(SNR) £ Q. (23) 

Thus, differential evolution [25] can then be employed to find 
the protograph ensemble with the lowest threshold. We next 
instantiate the iterative decoding threshhold for BEC surrogates 
and biAWGN surrogates. 

Definition 1 (Iterative decoding threshold under BEC surro¬ 
gate). The iterative decoding threshold of a protograph under 
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Figure 6 . Densities of L-values for 8 -ASK with shaped inputs after 100 iterations at SNR = 8.1 dB. 


the BEC surrogate is the minimum SNR on the original 
channel such that the resulting erasure probability vector 
£ = (ey( 1 ), €t( 2)) • • •) e T(N)) °f the BEC surrogates belongs 
to the protograph convergence region C e . 

Definition 2 (Iterative decoding threshold under the biAWGN 
surrogate). The iterative decoding threshold of a protograph 
under biAWGN surrogate is the minimum SNR on the original 
channel such that the resulting LLR standard deviation vector 
cTch = (^ ch ,T(i),c r ch , T( 2 ),...,cr ch , r(Ar) ) of the biAWGN sur¬ 
rogates belongs to the protograph convergence region C „ A . 

2 ) Shaped Input: For shaped input, the definition of the 
iterative decoding threshold is ambiguous. Since both the input 
distribution P B and the constellation scaling A influence both 
SNR and the bit-uncertainty set, there are infinitely many bit- 
uncertainty sets corresponding to the same value of SNR. We 
use an information theoretic argument to uniquely parametrize 
the bit-uncertainty set as a function of SNR. For each SNR, 
we choose P B and A that maximize the achievable rate, i.e., 
we let Pb(SNR) and A(SNR) be the solutions of 

max 1 H(B) — + Z) 1 . (24) 

A,P B : E[|Azb| 2 ]<SNR 1 I 

As in the uniform case, the bit uncertainty set is now 
parametrized by the SNR and we define the threshhold by 
(23). Definitions 1 and 2 now apply accordingly for shaped 
input. Note that for shaped input, we do not search along a 
constant rate but we search along the transmission rate. In 
Fig. 5, we display example search trajectories for uniform 
inputs and for shaped inputs for 8-ASK: In case of uniform 
inputs, the transmission rate is 2 bits/channel use for a rate 
(m — l)/m code independently of the SNR (dotted curve), 
whereas the transmission rate is given by H(B) — 1 in case of 
shaped input (dashed curve). In the following, we must keep 
in mind that for shaped input, changing the SNR also changes 
the transmission rate. 

V. Simulation Results of the Optimized Codes 

We choose the basematrix A from the set 
{0,1 ,...,S} MxN . The parameters M and N are chosen 


such that the desired code rate c = ( N — M)/N is met 
and the number of different degrees per bit-level D (see 
IV-B) is guaranteed. The maximum number of parallel edges 
S is crucial for both the performance of the code and the 
success of the optimization procedure. The product M ■ S 
describes the maximum variable node degree in the final code 
and thereby determines the size of the optimization search 
space. The optimization of the basematrices is performed as 
described in IV-C. 

A. Properness of BEC and biAWGN Surrogates 

We designed codes using both BEC and biAWGN surro¬ 
gates. Our simulations results reveal that BEC surrogates are 
not proper, i.e., the performance of the obtained codes on the 
original channel can be improved further by changing the bit- 
mapper. For the codes we obtained by using biAWGN surro¬ 
gates, the resulting bit-mapper are optimal, i.e., changing the 
bit-mapping degrades the performance on the original channel. 
Thus, for the considered scenario, the biAWGN channel is a 
proper surrogate and can be used for code design. 

This observation is also verified by comparing the densities 
of the L-values of both the original bit-channels and the 
biAWGN surrogates after decoding as illustrated in Fig. 6. 
The plot shows the densities Pl,|B;(^|0) an d P\_ |b. W0) (see 
IV-B2) for 8-ASK with shaped input at SNR = 8.1 dB. We 
observe a close match of the two densities. 

To further illustrate the close matching of the biAWGN 
surrogates with the original channel, we also plot simulation 
results for the biAWGN surrogates, see the curves with triangle 
markers in Fig. 9, Fig. 10 For each SNR, we calculate the 
corresponding set of surrogate channel parameters <r^ SNR) 
and simulate transmission over the corresponding uniform 
input symmetric-output biAWGN channels. 

B. Code Design for Uniform Input 

In order to compare our code design approach to the 
setting of Zhang and Kschischang [10], we design codes 
of rates 1/2 and 3/4 for 4-ASK constellations with uniform 
inputs. The optimized protographs are depicted in Table I. 
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Table I. OPTIMIZED PROTOGRAPHS AND DESIGN PARAMETERS FOR UNIFORM INPUT. 


Constellation (code rate) Basematrix 


P-EXIT Gap of Threshold to 

Decoding -Rbmd (5) 

Threshold 


/211214\ 

4-ASK (1/2) A= 1 1 1 2 2 5 5.57dB 0.28dB 

\1 00 1 06y 

B 2 Bi 

4-ASK (3/4) Aes (ll 226622 ) 9.57 dB 0.26dB 

b 2 b 1 


Table II. OPTIMIZED PROTOGRAPHS AND DESIGN PARAMETERS FOR SHAPED INPUT. 


Constellation (code rate) Basematrix 


8-ASK (2/3) 


/ 1 1 1 1 1 6 \ 

( 2 2 1 1 2 6 ) 

B 2 B 3 Bi 


P-EXIT Gap of Threshold to 

Decoding Ilog 2 (l + SNR) 

Threshold 


7.74 dB 0.39 dB 


64-ASK (5/6) 


A = 


222122622066\ 

111211610266; 



25.52 dB 


0.35 dB 


We discuss the rate 1/2 code in more detail. The coded bits 
are transmitted over two different bit-channels. For each bit- 
channel, we allow D = 3 (possibly different) variable degrees 
so that always 3 variable nodes in the protograph are assigned 
to the same surrogate channel: cr c h,i = er c h,2 = <r c h,3 and 
rich ,4 = cr c h, 5 = cr C h,6- Once the optimized basematrices have 
been found, we construct quasi-cyclic parity-check matrices 
with blocklengths n = 16200. We simulate the constructed 
codes using 100 decoding iterations. The bit error rates (BER) 
and frame error rates (FER) in Fig. 7 and 8 show that the 
finite length performance of our codes is equal to or slightly 
better than the codes in [10], As a reference, we also plot 
the Shannon limit for BMD and a 4-ASK constellation for the 
corresponding code rates 1/2 and 3/4. The decoding thresholds 
given in Table I are obtained by PEXIT analysis for the 
surrogate channels. 

C. Code Design for Shaped Input 

In the following we design codes of rate 2/3 for 8-ASK and 
rate 5/6 for 64-ASK to show that our design technique scales 
with the number of bit-levels. The optimized base matrices 
can be found in Table II. In order to limit the search space 
for the optimization, we choose D = 2 and allow a maximum 
of S = 6 parallel edges. For comparison, we use the off-the- 
shelf DVB-S2 LDPC codes used in [21] and use the optimized 
bit-mappers for DVB-S2 codes listed in [21, Table V]. 

As the shaped input causes a non-constant transmission rate 
as shown in IV-C2, the bit and frame error ratios are plotted 
over the gap to continuous AWGN capacity, i.e., 

ASNR| dB = SNR(FER/BER,C)| dB - 101og 10 (2 2 i? - l) , 

(25) 

where SNR(FER/BER,C)| dB denotes the SNR in dB that is 
necessary to achieve the desired FER/BER with code C and 



Figure 7. 4-ASK, uniform inputs: Performance of a rate 1/2 optimized 
protograph code compared to the rate 1/2 MET code in [10]. 


the rate R equals the transmission rate H (B) —1 at this specific 
SNR. We use the same input distributions for our optimized 
codes and for the DVB-S2 codes. 

As for the uniform case, the optimized protograph base- 
matrices have been lifted yielding quasi-cyclic parity-check 
matrices of blocklengths n = 64800. 100 decoding iterations 
have been performed. The simulation results are depicted in 
Fig. 9 and 10: At a target FER of HP 3 , the rate 2/3 code 
has a gap of 0.63dB to the AWGN capacity \ log 2 (l + SNR) 
at a spectral efficiency of 1.38 bits/channel use. For the 5/6 


























Figure 8. 4-ASK, uniform inputs: Performance of a rate 3/4 optimized 
protograph code compared to the rate 3/4 MET code in [ 0]. 



Figure 9. 8-ASK, shaped inputs: Performance comparison of optimized rate 2/3 
code and bit-mapping optimized DVB-S2 code. The horizontal axis denotes 
the gap to AWGN capacity as explained in (25). 

rate code, we observe a gap of 0.69 dB compared to AWGN 
capacity at a target FER of 10“ 3 and a spectral efficiency of 
4.25 bits/channel use. Hence, the bit-mapping optimized DVB- 
82 codes are outperformed by 0.1 dB (2/3 code) and 0.25 dB 
(5/6 code). 

VI. Conclusion 

We proposed a protograph-based LDPC code design ap¬ 
proach for bandwidth-efficient coded modulation that is suit¬ 
able both for uniform and shaped inputs. The different bit- 



Figure 10. 64-ASK, shaped inputs: Performance comparison of rate 5/6 code 
and bit-mapping optimized DVB-S2 code. The horizontal axis denotes the gap 
to AWGN capacity as explained in (25). 

channels are replaced by biAWGN surrogates so that PEXIT 
and differential evolution give ensembles with good decoding 
thresholds. The performance of the new codes for uniform 
inputs are as good as the best codes in literature. For shaped 
inputs, the new codes operate within 0.7dB to ilog 2 (l + 
SNR). Future research should investigate the influence of 
the surrogates on the code performance by employing a 
full-fledged density evolution for protographs. Furthermore, 
precoded protographs should be considered to improve the 
threshold without increasing the variable node degrees. 
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